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Abstract 

We consider the most general three point function for gravitational waves produced during a 
period of exactly de Sitter expansion. The de Sitter isometries constrain the possible shapes to 
only three: two preserving parity and one violating parity. These isometries imply that these 
correlation functions should be conformal invariant. One of the shapes is produced by the ordinary 
gravity action. The other shape is produced by a higher derivative correction and could be as large 
as the gravity contribution. The parity violating shape does not contribute to the bispectrum 
[l], [2I, even though it is present in the wavefunction. We also introduce a spinor helicity formalism 
to describe de Sitter gravitational waves with circular polarization. 

These results also apply to correlation functions in Anti-de Sitter space. They also describe the 
general form of stress tensor correlation functions, in momentum space, in a three dimensional 
conformal field theory. Here all three shapes can arise, including the parity violating one. 
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1. INTRODUCTION 



Recently there has been some effort in understanding the non-gaussian corrections to 
primordial fluctuations generated during inflation. The simplest correction is a contribution 
to the three point functions of scalar and tensor fluctuations. For scalar fluctuations there is 
a classification of the possible shapes for the three point function that appear to the leading 



orders in the derivative expansion for the scalar field [S-Ig]. 

In this paper we consider tensor fluctuations. We work in the de Sitter approximation 
and we argue that there are only three possible shapes for the three point function to all 
orders in the derivative expansion. Thus the de Sitter approximation allows us to consider 
arbitrarily high order corrections in the derivative expansion. The idea is simply that the 
three point function is constrained by the de Sitter isometrics. At late times, the interesting 
part of the wavefunction becomes time independent and the de Sitter isometrics act as the 
conformal group on the spatial boundary. We are familiar with the exact scale invariance, 
but, in addition, we also have conformal invariance. The conformal invariance fixes the 
three point functions almost uniquely. By "almost", we simply mean that there are three 
possible shapes allowed, two that preserve parity and one that violates parity. We compute 
explicitly these shapes and we show that they are the only ones consistent with the conformal 
symmetry. In particular, we analyze in detail the constraints from conformal invariance. In 
order to compute these three shapes it is enough to compute them for a simple Lagrangian 
that is general enough to produce them. The Einstein gravity Lagrangian produces one of 
these three shapes 3, ^- The other parity conserving shape can be obtained by adding a 
J term to the action, where W is the Weyl tensor. Finally, the parity violating shape 
can be obtained by adding J W'^W, where W is the Weyl tensor with two indices contracted 
with an e tensor. The fact that the gravitational wave expectation value is determined by the 
symmetries is intimately connected with the following fact: in four dimensional flat space 
there are also three possible three point graviton scattering amplitudes [9i\ ^. Though the 
parity violating shape is contained in the wavefunction of the universe (or in related AdS 
partition functions), it does not arise for expectation values [l, 2|.^. Thus for gravitational 



^ In flat space, one has to complexify the momenta to have nonvanishing three point amplitudes. In de 

Sitter, they are the natural observables. 
^ The previous version of this paper incorrectly stated that the parity violating shape contributed to de 
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wave correlators in ^5' we only have two possible shapes, both parity conserving. 

We show that, under general principles, the higher derivative corrections can be as large 
as the term that comes from the Einstein term, though still very small compared to the 
two point function. In fact, we expect that the ratio of (777)/(77)^ is of order one for 
the ordinary gravity case, and can be as big as one for the other shape. When it becomes 
one for the other shape it means that the derivative expansion is breaking down. This 
happens when the scale controlling the higher derivative corrections becomes close to the 
Hubble scale. For example, the string scale can get close to the Hubble scale. Even though 
ordinary Einstein gravity is breaking down, we can still compute this three point function 
from symmetry considerations, indicating the power of the symmetry based approach for 
the three point function. This gravity three point function appears to be outside the reach 
of the experiments occurring in the near future. We find it interesting that by measuring 
the gravitational wave three point function we can directly assess the size of the higher 
derivative corrections in the gravity sector of the theory. Of course, there are models of 

nn 

inflation where higher derivatives are important in the scalar sector p^, and in that case 
too, the non-gaussian corrections are a direct way to test those models jsl, Q, 13 1. The 
simplicity of the results we find here is no longer present when we go from de Sitter to an 
inflationary background. However, the results we find are still the leading approximation in 
the slow roll expansion. Once we are away from the de Sitter approximation, one can still 



study the higher derivative corrections in a systematic fashion as explained in 

Our results have also a "dual" use. The computation of the three point function for grav- 
itational waves is mathematically equivalent to the computation of the three point function 
of the stress tensor in a three dimensional conformal field theory. This is most clear when 
we consider the wavefunction of the universe as a function of the metric, expanded around 
de Sitter space at late times jsl. This is a simple consequence of the symmetries, we are not 
invoking any duality here, but making a simple statement ^. From this point of view it is 
clear why conformal symmetry restricts the answer. If one were dealing with scalar opera- 



Sitter expectation values. The fact that there is no parity violation in the gravity wave bispcctrum was 
shown in [l[ Q • 

^ Of course, this statement is consistent with the idea that such a wavefunction can be computed in terms 
of a dual field theory. Here we are not making any assumption about the existence of a dual field theory. 



Discussions of a possible dual theory in the de Sitter context can be found in 15|, |16[ . 
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tors, there would be only one possible three point function. For the stress tensor, we have 
three possibilities, two parity conserving and one parity violating. The parity conserving 
three point functions were computed in 17|. Here we present these three point functions 
in momentum space. Momentum space is convenient to take into account the conservation 
laws, since one can easily focus on the transverse components of the stress tensor. However, 
the constraints from special conformal symmetry are a little cumbersome, but manageable. 
We derived the explicit form of the special conformal generators in momentum space and 
we checked that the correlators we computed are the only solutions. In fact, we found it 
convenient to introduce a spinor helicity formalism, which is similar to the one used in flat 
four dimensional space. This formalism simplifies the algebra involving the spin indices and 
it is a convenient way to describe gravitational waves in de Sitter, or stress tensor correlators 
in a three dimensional conformal field theory. In Fourier space the stress tensor has a three 
momentum k, whose square is non-zero. The longitudinal components are determined by 
the Ward identities. So the non-trivial information is in the transverse, traceless compo- 
nents. The transverse space is two dimensional and we can classify the transverse indices in 
terms of their helicity. Thus we have two operators with definite helicity, T^ik). In terms of 
gravitational waves, we are considering gravitational waves that have circular polarization. 
These can be described in a convenient way by defining two spinors A and A, such that 
= (k^ \k\)°'^. In other words, we form a null four vector, and we proceed as in the four 
dimensional case. We only have S'0(3) symmetry, rather than 5*0(1,3), which allows us to 
mix dotted and undotted indices. We can then write the polarization vectors as oc cr*j^A"A* 
(no bar), etc. This leads to simpler expressions for the three point correlation functions of 
the stress tensor in momentum space. We have expressed the special conformal generator in 
terms of these variables. One interesting aspect is that this formalism makes the three point 
function completely algebraic (up to the delta function for momentum conservation). As 
such, it might be a useful starting point for computing higher point functions in a recursive 
fashion, both in dS and AdS. This Fourier representation might also help in the construction 
of conformal blocks. The connection between bulk symmetries and the conformal symmetry 



18|-l22|. 



on the boundary was discussed in the inflationary context in 

The idea of using conformal symmetry to constrain cosmological correlators was also 
discussed in {23!], which appeared while this paper was in preparation. Though the point of 
view is similar, some of the details differ. In that paper, scalar fluctuations were considered. 
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However, scalar fluctuations, and their three point function, crucially depend on departures 
from conformal symmetry. It is likely that a systematic treatment of such a breaking could 
lead also to constraints, specially at leading order in slow roll. On the other hand, the 
gravitational wave case, which is discussed here, directly gives us the leading term in the 
slow roll expansion. 

The paper is organized as follows. In section |2] we perform the computation of the most 
general three point function from a bulk perspective. We also discuss the possible size of the 
higher derivative corrections. In section [3] we review the spinor helicity formalism in 4D flat 
spacetime, and propose a similar formalism that is useful for describing correlators of CFTs 
and expectation values in dS and AdS. We then write the previously computed three point 
functions using these variables. In section H] we review the idea of viewing the wavefunction 
of the universe in terms of objects that have the same symmetries as correlators of stress 
tensors in CFT. We also emphasize how conformal symmetry constrains the possible shapes 
of the three point function. In section |5] we explicitly compute the three point function for 
the stress tensor for free field theories in 3D, and show that, up to contact terms, they have 
the same shapes as the ones that do not violate parity, computed from the bulk perspective. 
The appendices contain various technical points and side comments. 

Note added: We have revised this paper correcting statements regarding the parity 
violating terms in the bispectrum, in light of [l, 2|. 



2. DIRECT COMPUTATION OF GENERAL THREE POINT FUNCTIONS 

In this section we compute the three point function for gravitational waves in de Sitter 
space. We do the computation in a fairly straightforward fashion. In the next section we 
will discuss in more detail the symmetries of the problem and the constraints on the three 
point function. 



2.1. Setup and review of the computation of the gravitational wave spectrum 



The gravitational wave spectrum in single- field slow roll inflation was derived in [24 1. 
Here we will compute the non-gaussian corrections to that result. As we discussed above, 
we will do all our computations in the de Sitter approximation. Namely, we assume that 
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we have a cosmological constant term so that the background spacetime is de Sitter. There 
is no inflaton or scalar perturbation in this context. This approximation correctly gives the 
leading terms in the slow roll expansion. We leave a more complete analysis to the future. 
It is convenient to write the metric in the ADM form 

^00 = -iV' + 9ijN'N^, go^ = gijN^ , Qi, = e^""' exp{j)ij (2.1) 

Where H is Bubble's constant. and Ni are Lagrange multipliers (their equations of motion 
will not be dynamical), and parametrizes gravitational degrees of freedom. The action 
can be expressed as 

S=^J V^{R - QH') = ^J^s {NR^'^ - 6NH' + N-\E,,E^^ - {E\f)) (2.2) 

Where Eij = ^{gij — ViNj — WjNi) and we define Mpf = StcGj^. We fix the gauge by 
imposing that gravity fluctuations are transverse traceless, 'ju = and di'jij = 0. Up to 
third order in the action, we only need to compute the first order values of the Lagrange 



multipliers and A^j js]. By our gauge choice, these are A^ = 1 and Aj = as there cannot 
be a first order dependence on the gravity fluctuations. Expanding the action up to second 
order in perturbations we find 

S2 = ^J [e^'^'^^r^ - e"%%A7^,) (2.3) 

We can expand the gravitational waves in terms of polarization tensors and a suitable choice 
of solutions of the classical equations of motion. If we write 

E <,e^-lc^{t)al, + h.c. (2.4) 

the gauge fixing conditions imply that the polarization tensors are traceless, en = 0, and 
transverse kieij = 0. The helicities can be normalized by efje*^ = 46^^ . The equations of 
motion are then given by 

2 e"-^* 
=i:i{v) - -Iciiv) + k'lciiv) , V = —rr (2-5) 

Tj £1 

where we introduced conformal time, r]. We take the classical solutions to be those that 
correspond to the Bunch-Davies vacuum [25^, so "ydiv) = "^^^^^^^^^ ~ ikrf). Here we have 
denoted by A; = \k\ the absolute value of the 3-momentum of the wave. We are interested in 
the late time contribution to the two-point function, so we take the limit where 77 —t- 0. After 



Fourier transforming the late-time dependence of the two-point function and contracting it 
with polarization tensors of same helicities we find: 

In the inflationary context (with a scalar field), higher derivative terms could give rise 
to a parity breaking contribution to the two point function. This arises from terms in the 
effective action of the form j f{(l))WW QjH]- This parity breaking term leads to a different 
amplitude for positive and negative helicity gravitational waves, leading to a net circular 
polarization for gravitational waves. If / is constant this term is a total derivative and it 
does not contribute. Thus, in de Sitter there is no contribution from this term. In other 
words, the parity breaking contribution is proportional to the time derivative of /. Some 
authors have claimed that one can get such parity breaking terms even in pure de Sitter 
28|. However, such a contribution would break CPT. Naively, we would expect that CPT 
is spontaneously broken because of the expansion of the universe. However, in de Sitter we 
can go to the static patch coordinates where the metric is static. For such an observer we 
expect CPT to be a symmetry. A different value of the left versus right circular polarization 
for gravitational waves would then violate CPT. 

In the AdS case, or in a general CFT, there can be parity violating contact terms in the 
two point function^. This is discussed in more detail in appendix O 



2.2. Three point amplitudes in flat space 

In order to motivate the form of the four dimensional action that we will consider, let 
us discuss some aspects of the scattering of three gravitational waves in fiat space. This is 
relevant for our problem since at short distances the spacetime becomes close to fiat space. 

In fiat space we can consider the on shell scattering amplitude between three gravita- 
tional waves. Due to the momentum conservation condition we cannot form any non-zero 
Mandelstam invariant from the three momenta. Thus, all the possible forms for the ampli- 
tude are exhausted by listing all the possible ways of contracting the polarization tensors 

^ A contact term is a contribution proportional to a delta function of the operator positions. 
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of the gravitational waves and their momenta, js]^. There are only two possible ways of 
doing this, in a parity conserving manner. One corresponds to the amplitude that comes 
from the Einstein action. The other corresponds to the amplitude we would get from a 
term in the action that has the form W^, where W is the Weyl tensor. In addition, we can 
write down a parity violating amplitude that comes from a term of the form WW"^, where 
Wabcd = ^abefW^^^^. Thcsc tcrms involving the Weyl tensor are expected to arise from higher 
derivative corrections in a generic gravity theory. By using field redefinitions, any other 
higher derivative interaction can be written in such a way that it does not contribute to the 
three point amplitude. 

By analogy, in our de Sitter computation we will consider only the following terms in the 
gravity action 



Seff = I d X 



-g [-f- + R)+ (a W'''',M'''mnW"'\b) 



^ (2.7) 

Here A is a scale that sets the value of the higher derivative corrections. We will discuss 
its possible values later. This form of the action is enough for generating the most general 
gravity three point function that is consistent with de Sitter invariance. This will be shown 
in more detail in section Hj by using the action of the special conformal generators. For the 
time being we can accept it in analogy to the fiat space result. Instead of the Weyl tensor in 
(12. 7p we could have used the Riemann tensor. The disadvantage would be that the term 
would not have vanished in a pure de Sitter background and it would also have contributed 
to the two point function. However, these extra contributions are trivial and can be removed 
by field redefinitions. So it is convenient to consider just the term. 



2.3. Three-point function calculations 

In this subsection we compute the three point functions that emerge from the action in 
(12. 7p . First we compute the three point function coming from the Einstein term, and then 
the one from the term. 

^ In flat space, the three point amphtude is non-trivial only after analytically continuing to complex values 
of the momentum. 
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2.4. Three point function from the Einstein term 

This was done in [3, O]^- For completeness, we review the calculation and give some 
further details. To cubic order we can set = 1, A^j = in (12. ip . Then the only cubic 
contribution from (12. 2p comes from the term involving the curvature of the three dimensional 
slices, R^^\ Let us see more explicitly why this is the case. On the three dimensional slices 
we define g = e'^^^g, with cjij = {e'^)ij- AH indices will be raised and lowered with g. The 
action has the form 

= 1 y dtd^x [e^*^(3) ^ e-'^'iEijE'^ - {E\f)\ (2.8) 

Now we prove that the second term does not contribute any third order term to the action. 
To second order in 7 we have = (7 + |[7, 7])jj- Then we find 

E,jE^^ - {E\Y = 7.,7^, + 0(7') (2-9) 

which does not have any third order term. Thus, the third order action is proportional to 
the curvature of the three-metric. This is then integrated over time, with the appropriate 
prefactor in (12. 8p . Of course, if we were doing the computation of the flat space three point 
amplitude, we could also use a similar argument. The only difference would be the absence 
of the e^* factor in the action (12.81) . Thus, the algebra involving the contraction of the 
polarization tensors and the momenta is the same as the one we would do in fiat space (in 
a gauge where the polarization tensors are zero in the time direction). Thus the de Sitter 
answer is proportional to the flat space result, multiplied by a function of \ki\ only, which 
comes from the fact that the time dependent part of the wavefunctions is different. 

We want to calculate the tree level three-point function that arises from this third- 
order action. In order to do that, we use the in-in formalism. The general prescrip- 
tion is that any correlator is given by the time evolution from the "in" vacuum up 
to the operator insertion and then time evolved backwards, to the "in" vacuum again, 
{0{t)) = (in |f e"*-^^'"*(*')'^*'0(t)Te*-^^™*(*')'^*'| in). We are only interested in the late-time 
limit of the expectation value. We find 

/+00 
dt' [H,nt{t'), 7'H^i,+oo)7'^(x2,+oo)7^«(x3,+oo)] 
■00 

(2.10) 



1 



In m the AdS case was considered. 



We write the gravitational waves in terms of oscillators as in (12. 4p . We calculate correlators 
for gravitons of specific helicities and 3- momenta. Note that, because there are no time 
derivatives in the interaction Lagrangian, then it follows that Hf^^ = —Lifnt- Once we put in 
the wavef unctions, the time integral that we need to compute is of the form /m[J°^ drf^il — 
ikir])(l — ik2r]){l — ik3r])e^^'''^~^''^~^''^^^] (in conformal time). Two aspects of the calculation are 
emphasized here. One is that we need to rotate the contour to damp the exponential factor 
at early times, which physically corresponds to finding the vacuum of the interacting theory 
8|, as is done in the analogous flat space computation. Another aspect is that, around zero, 
the primitive is of the form _e^l^ll^il^;£^ _ _i _ ^(/^..^ _|_ /^-^ + k^) + 0(e), where e is our 
late-time cutoff. This divergent contribution is real and it drops out from the imaginary 
part. We get 

i^tHWn - (2.)3.3 ^ + k3) (^)' 

mk]eiy,,el, - 2el^{kfel){klei^) + cyclic] x (2.11) 
fi- _L <i- ^1^2 + kiks + k2k^ kik2k3 \ 

r+ ^+ ' kr + k2 + k, JMThTW) 

The second line is the one that is the same as in the fiat space amplitude. The third 
line comes from the details of the time integral. Below we will see how this form for the 
expectation value is determined by the de Sitter isometrics, or the conformal symmetry. 

2.5. Three point amplitude from W'^ in flat space 

Let us calculate the following term in fiat space, to which we will refer as W^: 
W"-^ ^sW^^ apW'^^ ai3- We can write the following first order expressions for the components 
of the Weyl tensor 





2%- 


W'\k = 


2 i.iki,j - ikj,i) 






W'hi = 


2 {-^ikljl + ^illik + ^jklil - ^jllik) 



(2.12) 
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where we used that 7 is an on shell gravitational wave. i.e. 7 obeys the flat space equations 
of motion. We also used that — di^ij — 0, Ni — 0, N — 1. We can then write 

(2.13) 

Evaluating these terms leads us to 

s^^^ = y A"^ [27ij7jfc7fei + 'iii3iki,iiki,3 + ^inik,nji,k - Q^ij^ik,nki,j] (2.i4) 

Plugging 7y = eljC^'^^'^ + e^^e'^^^ _|_ ^s^.gifcs :r where k ■ x — k^Xi — kt, we get the following 
expression for the vertex due to the term: 

Vw^jiat = 6 kik2k3 [ki k'^kiel^elieli + cychc - 

(2.15) 

(A;i + A;2 + h){elklel,k{e% + cychc) - 2 kik2hele%el,] 
By choosing a suitable basis for the polarization tensors, one can show that this agrees with 
the gauge invariant covariant expression V^ijiat = Qk'lk'(e^^k2k2e^^k^^k^e^^. 

2.6. Three point function from in dS 

The straightforward way of performing the computation would be to insert now the 
expressions for the wavefunctions in the term in de Sitter space, etc. There is a simple 
observation that allows us to perform the de Sitter computation. First we observe that the 
Weyl tensor is designed so that it transforms in a simple way under overall Weyl rescaling 
of the metric. Thus the Weyl tensor for the metric in conformal time is simply given by 
Wfj_,^Sa{g) = -jjT^W ^iu5u{.g = ^^)- Note also that, for this reason, the Weyl tensor vanishes in 
the pure de Sitter background. Thus, we only need to evaluate the Weyl tensor at linearized 
order'^. For on shell wavefunctions 7 = (1 — ikr})e^^'^'^'^^-^ we can show that 

where VF-^'"* is the expression for the flat space Weyl tensor that we computed in the previous 
section, computed at linearized order for a plane wave around flat space. Thus, when we 

insert these expressions in the action we have 

fw^^ f dr]d^x{kxk2hrf){H''rf){W^^'''f (2.17) 
J J —00 



Note that the term does not contribute to the two point function. 
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The whole algebra involving polarization tensors and momenta is exactly the same as in 
flat space. The only difference is the time integral, which now involves a factor of the form 
J drjTfe^^'^ oc where we have defined E = ki + k2 + k^, and we rotated the contour 

appropriately. Putting all this together, we get the following result for the three-point 
function due to the term 



X 



H y fH^ (-30) (2.18) 



MpiJ \AJ {ki + k2 + k3)%kik2k3 

where Vwsjiat was introduced in fl2.15p . There are also factors of 1/kf that were included 
to get this result. The parity violating piece will be discussed after we introduce spinor 
variables, because they will make the calculation much simpler. 



2.7. Estimating the size of the corrections 

Let us write the effective action in the schematic form 

M2 ' 



S 



'PI 



J[^R-6H'^ + L' J W 



+ ■■■ (2.19) 



2 

where the dots denote other terms that do not contribute to the three point function. Here 
L is a constant of dimensions of length. We have pulled out an overall power of Mpi for 
convenience. The gravitational wave expectation values coming from this Lagrangian have 
the following orders of magnitude 

H'^ H'^ 
(77) - ^ , (777)k = ^ (777)h^3 = j^iLH)' (2.20) 

Thus the ratio between the two types of non-gaussian corrections is 

(777) H/3 .4 0-4 /oon 
^ L H [2.21) 

(777) 

We know that H^/Mpi is small. This parameter controls the size of the fluctuations. In 
the AdS context, we know that when the right hand side in fl2.2ip becomes of order one we 



have causality problems 29l-l32| . We expect that the same is true in dS, but we have not 
computed the precise value of the numerical coefficient where such causality violation would 
occur. So we expect that 

HL < 1 (2.22) 
12 



In a string theory context we expect L to be of the order of the string scale, or the Kaluza 
Klein scale. Thus the four dimensional gravity description is appropriate when HL ^ 1. 
In fact, in string theory we expect important corrections when Hlg ~ 1- In that case, the 
string length is comparable to the Hubble scale and we expect to have important stringy 
corrections to the gravity expansion. Note that in the string theory context we can still 
have H'^/Mpi ~ g"^ being quite small. So we see that there are scenarios where the higher 
derivative corrections are as important as the Einstein contribution, while we still have a 
small two point function, or small expansion parameter H'^/Mpi. In general, in such a 
situation we would not have any good argument for neglecting higher curvature corrections, 
beyond the term. However, in the particular case of the three point function, we can 
just consider these two terms and that is enough, since these two terms (the Einstein term 
and the term) are enough to parametrize all the possible three point functions consistent 
with de Sitter invariance. If we define an fNL-gravUy = then we find that the 

Einstein gravity contribution of fNL-gravUy is of order one. This is in contrast to the f^L 
for scalar fluctuations which, for the simplest models, is suppressed by an extra slow roll 
factor^. 

In an inflationary situation we know that the fact that the fluctuations are small is 
an indication that the theory was weakly coupled when the fluctuations were generated. 
However, it could also be that the stringy corrections, or higher derivative corrections were 
sizable. In that case, we see that the gravitational wave three point function (or bispectrum) 
gives a direct measure of the size of higher derivative corrections. Other ways of trying to 



see these corrections, discussed in SSj , involves a full reconstruction of the potential, etc. In 
an inflationary context terms involving the scalar field and its time variation could give rise 
to new shapes for the three point function since conformal symmetry would then be broken. 
However, one expects such terms to be suppressed by slow roll factors relative to the ones we 
have considered here. However a model specific analysis is necessary to see whether terms 
that contain slow roll factors, but less powers of LH dominate over the ones we discussed. For 
example, a term of the form MpiL'^f{(f))W'^ is generically present in the effective action jsj. 
Such a term could give a correction of the order {m) fw"^ / {ill) r ~ ^f^HL^. Here ej is 



Note that we have divided by the gravity two point function to define fNL-gravUy If we had divided by 
scalar correlators, we would have obtained a factor of e^. 
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a small quantity of the order of a slow roll parameter, involving the time derivatives of /. 
Whether this dominates or not relative to f l2.2ip depends on the details of the inflationary 
scenario. In most cases, one indeed expects it to dominate. It would be very interesting 
if (I2.2ip dominates because it is a direct signature of higher derivative corrections in the 
gravitational sector during inflation. 

Notice that the upper bound (12.221) is actually smaller than the naive expectation from 
the point of view of the validity of the effective theory. From that point of view we would 
simply demand that the correction due to at the de Sitter scale H should be smaller 
than one. This requires the weaker bound H^L'^ < This condition is certainly too lax 
in the AdS context, where one can argue for the more restrictive condition fl2.22p . 

In summary, we can make the higher derivative contribution to the gravity three point 
function of the same order as the Einstein Gravity contribution. Any of these two terms 
are, of course, fairly small to begin with. 



3. SPINOR HELICITY VARIABLES FOR DE SITTER COMPUTATIONS 



In this section we introduce a technical tool that simplifies the description of gravitons 
in de Sitter. The same technique works for anti-de Sitter and it can also be applied for 
conformal field theories, as we will explain later. 

The spinor helicity formalism is a convenient way to describe scattering amplitudes of 



bur dimensions. We review the basic ideas here. For a more 



massless particles with spin in 
detailed description, see Isf 



36|. In four dimensions the Lorentz group is SO (1,3) ~ 
SL{2) X SL{2). A vector such as A;^ can be viewed as having two SL{2) indices, k"-^. The 
new indices run over two values. A 4-momentum that obeys the mass shell condition, 
= can be represented as a product of two (bosonic) spinors k"-^ = A"A^. Note that if we 
rescale A — )■ wX and A — i- ;^A we get the same four vector. We shall call this the "helicity" 
transformation. Similarly, the polarization vector of a spin one particle with negative 
helicity can be represented as 

r"^ = 7^ (3.1) 

where we used the SL{2) invariant contraction of indices (A,yu) = eabX°'fJ'^, where eab is the 
SL(2) invariant epsilon tensor. We have a similar tensor e^^ to contract the dotted indices. 
We cannot contract an undotted index with a dotted index. Note that this polarization vec- 
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tor (13.11) is not invariant under the helicity transformation. In fact, we can assign it a definite 
helicity weight, which we call minus one. This polarization tensor (13. ip is independent of 
the choice of ft. More precisely, different choices of ft correspond to gauge transformations 
on the external particles. For negative helicity we exchange A, f/ <-> A, r/ in (13. ip . For the 
graviton we can write the polarization tensor as a "square" of that of the vector 

..a,.b\a\b \a \b j-.d j-.b 

C+abab _ ^ ^-abab _ ^ /o 

^ " ' ^ " {Xn? ^ ^ ^ 

The product of two four vectors can be written as k.k' = — 2(A, A') (A, A'). 

Now let us turn to our problem. We are interested in computing properties of gravitational 
waves at late time. We still have the three momentum k. This is not null. However, we can 
just define a null four momentum {\k\,k). This is just a definition. We can now introduce A 
and A as we have done above for the flat space case. In other words, given a three momentum 
k we define A, A via 

{\k\,ky^ = {\k\ + k.a^'') = A"A^ (3.3) 

In the de Sitter problem we do not have full SL{2) x SL{2) symmetry. We only have one 
SL{2) symmetry which corresponds to the 50(3) rotation group in three dimensions. This 
group is diagonally embedded into the SL{2) x SL{2) group we discussed above. In other 
words, as we perform a spatial rotation we change both the a and a indices in the same way. 
This means that we now have one more invariant tensor, e^^ which allows us to contract the 
dotted with the undotted indices. For example, out of A'^ and A^ we can construct (A, A) 
by contracting with e^^^. This is proportional to \k\. Thus, this contraction is equivalent 
to picking out the zero component of the null vector. When we construct the polarization 
tensors of gravitational waves, or of vectors, it is convenient to choose them so that their 
zero component vanishes. But, we have already seen that extracting the zero component 
involves contracting dotted and undotted indices. We can now then choose a special /2 in 
(13. ip which makes sure that the zero component vanishes. Namely, we choose Ji^ = \^ . This 
would not be allowed under the four dimensional rules, but it is perfectly fine in our context. 
In other words, we choose polarization vectors of the form 

^ (A, A) ' ^ (A, A) ^'-^^ 

Notice that the denominator is just what we were calling k = \k\. Note also that the zero 
component of ^ is zero, since this involves contracting the a and b indices. This gives a 
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vanishing result due to the antisymmetry of the inner product. In our case we have a delta 
function for momentum conservation due to translation invariance, but we do not have one 
for energy conservation. The delta function for momentum conservation can be written by 
contracting A^A^ with cr*"'' in order to get the spatial momentum. Alternatively we can 
say that AjAj oc e'*^. This is just saying that the fourvector has only a time component. 

For the graviton, we likewise take /i = A and /i = A in (13. 2p . With these choices we make 
sure that the polarization vector has zero time components and that it is transverse to the 
momentum. 

Everything we said here also applies for correlation function of the stress tensor in three 
dimensional field theories. If we have the stress tensor operator Tij{k) in Fourier space, 
we can then contract it with a polarization vector transverse to k constructed from A and 
A. In other words, we construct operators of the form T"*" = C,i~^j~Tij with C,~^ as in (13.41) . 
This formalism applies for any case where we have a four dimensional bulk and a three 
dimensional boundary, de Sitter, Anti-de Sitter, Hyperbolic space, Euclidean boundary, 
Lorentzian boundary, etc. The only difference between various cases are the reality condi- 
tions. For example, in the de Sitter case that we are discussing now, the reality condition is 
(A'^)* = e„,A^ 

In summary, we can use the spinor helicity formalism tyo describe gravitational waves in 
de Sitter, or any inflationary background. It is a convenient way to take into account the 
rotational symmetry of the problem. One can rewrite the expressions we had above in terms 
of these variables. 



3.1. Gravitational wave correlators in the spinor helicity variables 

Let us first note the form of the two point function. The only non-vanishing two point 

functions are the -|--|- and two point functions. This is dictated simply by angular 

momentum conservation along the direction of the momentum. Since the momenta of the 
two insertions are opposite to each other, their spins are also opposite and sum to zero as 
they should. The two point functions are then 

(7+7^) = + k')^-^ = S'ik + k')-l- (3.5) 
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where in the last formula we have used a particular expression for A' in terms of A. More 
precisely, if the momentum of one wave if k, with its associated A and A, then for k' = —k 
we can choose A' = A and A' = —A. Here we have used that the matrices cr*^^ are symmetric. 

In the first expression we can clearly see the helicity weights of the expression. For the 

one we get a similar expression. 

We can now consider the three point functions. The simplest to describe are the ones 

coming from the interaction. In fact, these contribute only to the + + + and 

correlators, but not to the + H — correlators. This is a feature which is also present in the 
flat space case. These non vanishing correlators can be rewritten as 



M = {27rY6' (ki + k2 + ka) 



MpiJ \A 

where the kn in the denominators can also be written in terms of brackets such as 
kn = —{n,n), if so desired. Note that, when rewritten in terms of the A„ and A„, the 
above expressions are just rational functions of the spinor helicity variables (up to the over- 
all momentum conservation delta function). One can check that indeed the + H — and 

h vertices vanish for the term, which is straightforward by using the expressions in 

appendix [Bl Note that this is not trivial because we do not have four-momentum conser- 
vation, only the three-momenta are conserved. The parity violating interaction W^W does 
not contribute to the de Sitter expectation values 

The Einstein term contributes to all polarization components 



~{k^k2 + k.ks + k2k3){k^ + k2 + k,) - k^k^k^] [(1,2)(2,3)(3, 1)]^ (3.7) 
{<lt%,)R = (2vr)^5^ (^Yl {-M^) \^kXk,)^ f^^^ + - ^^^^^^ " + ^^^^^^ + - 



kik2 + kik'i + /ca^s fcifc2fc3 

ki + k2 + k^- 



(1,2)= 



(1,3) (3, 2) 



(3.8) 



ki + k2 + k3 {ki + k2 + k3y, 

and similar expressions for h and . Note that the Einstein gravity contribution to 

+ + + or is non- vanishing. This is in contradistinction to what happens in fiat space, 
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where it does not contribute to the + + + or cases. This might seem surprising, given 

that we had said before that the polarization tensor contribution to the time integrand is 
the same as the flat space one. After doing the time integral, in flat space we get energy 
conservation, which we do not have here. This explains why we got a non-vanishing answer. 
In fact, the flat space amplitude is recovered from the above expressions by focusing on the 
coefficients of the double poles in E = ki + k2 + k^. The fact that (13. 7p does not have a 
double pole ensures that the flat space answer is zero for those polarizations^. Similarly, the 
flat space answers for are obtained by looking at the coefficient of the 6^^ order pole in 



The expressions (13. 7p can also be written in a form that shows explicitly the effect of 
changing the helicity of one particle: 



In the next section we will show that the forms of these results follow from demanding 
conformal symmetry. 

4. GRAVITATIONAL WAVE CORRELATION FUNCTION AND CONFORMAL 
SYMMETRY 

In this section we will show how the three point functions we discussed above are con- 
strained by conformal symmetry. 

4.1. Wavefunction of the universe point of view 

In order to express the constraints of conformal symmetry it is convenient to take the 
following point of view on the computation of the gravity expectation values. Instead of 

^ In comparing to the flat space result, there are also factors of {kik2kz) that come from the normalization 
of the wavefunction. 



E in f lXB]) . 




(3.9) 




(3.10) 



(3.11) 
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computing expectation values for the gravitational waves, we can compute the probability 
to observe a certain gravitational wave, or almost equivalently, the wavefunction \I^(7). The 
expectation values are given by simply taking |\E'(7)p and integrating over 7. This point 
of view is totally equivalent to the usual one, where one computes expectation values of 7. 
It is useful because it makes the connection to AdS very transparent^*^. It also makes the 
action of the symmetries more similar to the action of the symmetries in a conformal field 
theory. This is explained in more detail in [sl (see also 37|). 
One writes the wavefunction in the form: 

v^ = expfi [ d'xd'y{T%x)T^'iy))Yix)Y'iy) + 



1 ^ 

+-j d'xd^d'z{T^ix)T^'iy)r"iz))Yix)Y'iy)Y"iz) + 

The first term expresses the simple fact that the wavefunction is gaussian. From this point 
of view, the quantity (T'^(x)T*'(?/)) is just setting the variance of the gaussian. Namely, 
this is just a convenient name that we give to this variance. Similarly for the cubic term, 
which is responsible for the first non-gaussian correction, etc. Here we have ignored local 
terms that are purely imaginary and which drop out when we take the absolute value of the 
wavefunction. From this expression for the wavefunction one can derive the following forms 



for the two and three point functions js], to leading order in the loop expansion, 

/rpSirpS2rpS3\ j_ IrpSi rpS2 rpS^ \* 

l^s, S2 s,^ ^ KJ-k,-^k2-^k3l + -k^-^ -k2'^ -k,) ON 

\ /fci ik2 ik-3,1 n-(2(T*'T'*' .)) \ '^i 

So we see that it is easy to go from the description in terms of a wavefunction to the 
description in terms of expectation values of the metric. The complex conjugate arises from 
doing and we used that 7*(— A;)* = However, if the wavefunction contains terms 

that are pure phases, we can loose this information when we consider expectation values of 
the metric. Precisely this happens when we have the parity violating interaction J W^W. 
It contributes to a term that is a pure phase. 

Here \1/ is the usual Wheeler de Witt wavefunction of the universe, evaluated in pertur- 
bation theory. It is expressed in a particular gauge, because we have imposed the = 1, 



In fact, the perturbative de _Sitter computation is simply an analytic continuation of the perturbative 
Anti-de Sitter computation 
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Ni = conditions. The usual reparametrization constraints and Hamiltonian constraints 
boil down to some identities on the functions appearing in (14.11) . These identities are pre- 
cisely the Ward identities obeyed by the stress tensor in a three dimensional conformal field 
theory^^. In the AdS case, this is of course familiar from the AdS/CFT point of view. In 
the de Sitter case, it is also true since this wavefunction is a simple analytic continuation 
of the AdS one. It is an analytic continuation where the radius is changed by i times the 
radius. In any case, one can just derive directly these Ward identities from the constraints of 
General Relativity. These identities express the fact that the wavefunction is reparametriza- 
tion invariant. For the case that we have scalar operators (and corresponding scalar fields 
in dS) we get an identity of the form di{Tij{x) Ylk 0{xk)) = — Ylii ^'^(^ ~ ^')^x^ (IIa: 0{xk))- 
These are derived by starting with the reparametrization constraint, taking multiple deriva- 
tives with respect to the arguments of the wavefunction, and setting all fiuctuations to zero 
after taking the derivatives. There is also another identity coming from the Hamiltonian 
constraint. This involves the trace of T and it takes into account the dimension of the 
operator. Namely we have (Tjj(x) Y\kO{xk)) = ~ ^O^KIlfc ^(^fc))- From these 

two identities, we can derive equations for the correlation functions if we have a conformal 
Killing vector Alternatively, we can derive these equations simply by noticing that a con- 
formal reparametrization does not change the metric on the boundary, up to a rescaling (or 
a shift of time in the bulk). Thus, this leaves the wavefunction explicitly invariant, without 
even changing the metric, which is why we get equations on correlation functions for each 
isometry of the background space. From the general relativity point of view, this is just the 
statement that each isometry of the background leads to a constraint on the wavefunction. 
In our case, the operators are other insertions of the stress tensor. Thus, we can think of the 
coefficients (TT) and (TTT) appearing in (14. ip as correlation functions of "stress tensors" . 
We are not assuming the existence of a dual CFT, we are simply saying that these quantities 

Though the Ward identities are the same, some of the positivity constraints of ordinary CFT's are not 
obeyed. For example, the {TT) two point function is negative. Thus, if there is a dual CFT, it should 
have this unusual property. 
^'^ If is the conformal Killing vector, then we can multiply the diTij ■ ■ ■ equation by u^ , integrate over x, 
integrate by parts, use the conformal Killing vector equation = ^riij{d.v), use the Ta equation and 

obtain the equations '^g[v'^{xs)dj:i + ^■"(^°) As](nfc 0{xk)) — 0. These encode aU the equations obeyed 
by correlators that are a consequence of the de Sitter isometrics at late times. For the metric, or stress 
tensors, the equations contain more indices and we write them in detail below. 
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obey the same Ward identities as the ones for the stress tensor in a CFT. A more precise 
discussion of these identities can be found in appendix |E] and in section 14.21 

The isometrics of de Sitter translate into symmetries of the wavefunction. Some of these 
are simple, like translation invariance. A less trivial one is dilatation invariance, or scale 
invariance. This simply determines the overall scaling of the three point function in terms 
of the momentum. If we think of 7 as a dimensionless variable, then its Fourier components 
have dimension minus three. Thus the total dimension of any n point function is — Sra. The 
delta function of momentum conservation takes into account a —3, and the remainder is the 
overall degree of homogeneity in the momentum. For the two point function it is —3 and 
for the three point function it is —6. It is a simple matter to count powers of momenta in 
the expressions we have given in order to check that this is indeed the case. 

If instead we look at correlators of the stress tensor, then in position space, we have that 
the operator has dimension three, while in momentum space it has dimension zero. 

The constraints from special conformal transformation are harder to implement and we 
discuss them in the next section. The results of the following section are also valid in any 
three dimensional conformal field theory. The general form for the three point function in 
position space was given in Here we study the same problem in momentum space. The 
expressions we find seem a bit simpler to us than the ones in [l^, but the reader can judge 
by him or herself. 

4.2. Constraints from special conformal transformations 

de Sitter space is invariant under a full S'0(l,4) symmetry group. The metric ds^ = 
~'^^^t'^^ makes some of these isometrics manifest. In particular the scaling symmetry 
changes Xi — )■ axi and 1] — )■ arj. There are also three more isometrics that are given in 
infinitesimal form by 

X* + 6*(-r7^ + f.x) - 2x*(6.a;) , i] i] - 2r](b.x) (4.4) 

where b is infinitesimal. When 77 — 0, which is the future boundary of the space, the 
time rescaling acts in a simple way on the wavefunction. In addition we can drop the 77^ 
in the space part. The transformation then becomes what is called a "special conformal" 
transformation on the boundary parametrized by x. In this section we study the action of 
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these transformations in detail. 

Now we work with coordinates on the three dimensional boundary. A special conformal 
transformation is given by 



= x%' - 2x\x.h) (4.5) 

Y} ^ = = 2{x^V - x'V) - 25}(x.6) = 2Mh - 25]{x.h) (4.6) 

1 ^ 

J = det(l + S)i/3 ^ 1 + V S% = 1 - 2{x.h) (4.7) 

u=l 

where 6* is an infinitesimal parameter. The transformation law for a tensor is 

rK..<M = jsr^(^--^)r„„..(x) (4,8) 

where A is the conformal dimension. For a current or the stress tensor we have A = 2, 3 



respectively. These transformation laws describe the (infinitesimal) action of the de Sitter 
isometrics on the comoving coordinates at late time. 

Now, in order to compute the variation of a correlator, we are interested in its change as 
a function. This means that the transformed correlator, as a function of the new variables, 
x' should be the same as the old correlator as a function of x. Thus we can evaluate T'{x) 
(and not x'). Then we write x = x' — 5x. In that way we find that the change is 

n 

6T^,...„„ =A2{x.b)T^,...„„ -2Y,M'''^T^,...,,...^„ - DT^,...^„ 

1=1 

D = x^{b.d)-2{b.x){x.d) (4.9) 

The matrix M was defined in (14.61) . We now Fourier transform (14. 9 p . The terms that 
contain a single power of x are easy to transform. They are given simply by inserting factors 
of a; — > —idk- For the term involving a it is important that we first replace the x — > —idk 
and then we change the derivatives by factors of /c, dx — )■ —ik. Thus a term like 

x^di ^idlh = i{kidl + 29^.) (4.10) 
x\x.dx) -^i{,dkr{dk3kj)) = i[Adk: + kjdk^dk'] (4.11) 



22 



We then find that the special conformal generator (up to an overall i), now has the form 

n 

ST,,...,,Xk) = - (A - 3)2(6.9fc)T,,...,„(fc) + 2J2M^'^T^,■■■n■■■^nik) - 5T,,...,„(fc) 

1=1 

D={b.k)dl-2kjdk^{b.dk) (4.12) 

Where the (—3) in A — 3 comes from the commutators we had in fl4.10p . 

In momentum space any generator has an overall momentum conserving delta function 
S{Y2jki). It is possible to pull the momentum space operator through the delta function. 
One can show that all terms involving derivatives of the delta function vanish. This is argued 
in detail in appendix [Dl The final result is that we can simply act with the operator fl4.12p 
on the coefficient of the delta function. 

We would now like to express the action of the special conformal generator in terms of 
the spinor helicity variables. This problem is very similar to the one analyzed for amplitudes 
in [341 . There it was shown that the special conformal generator is given by 

b.d = h'a""^—— (4.13) 

The closure of the algebra implies that this simple form can only be consistent when it is 
applied to objects of scaling dimension minus one (in Fourier space). In our case, we will 
see that 04.131) differs from the special conformal generator only by terms proportional to 
the Ward identity for the corresponding tensor (the current or the stress tensor). This will 
be discussed in more detail below. 



4.3. Constraints of special conformal invariance on scalar operators 



The correlation function of three scalar operators is very simple in position space and it is 
given by a well known formula. In momentum space, it is hard to find an explicit expression 
because it is hard to do the Fourier transform. For the case of the three point function, the 
answer is a function of the In that case we can rewrite the special conformal generator 
as 



h.k 



(4.14) 
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We see that the case of A = 2 is particularly simple^^. So we consider a situation with 
three scalar operators of dimension A = 2. Invariance under special conformal transforma- 
tions then implies 

k\dlj + mj + kldlj = Q (4.15) 

where / is the Fourier transform of the correlator. Using momentum conservation we can 
conclude that all second derivatives should be equal, for any m ^ n: 

idl^-dl)fih,k2,h) = (4.16) 

For each pair of variables this looks like a two dimensional wave equation. Thus the general 
solution is given by /(fci, ^2, h) = 9ih + + h) + H^i - ^2 - h) + /(^2 -h- h) + m{k3 - 
ki — k2). To fix the form of these functions we look at the dilatation constraint: 

{kidk^ + k2dk^ + /c3(9fc3)/(/ci, k2, k'i) = c (4.17) 

In principle, c = 0. We would be tempted to conclude that this implies that each of the 
functions in / should be scaling invariant. This would leave only a constant solution. One 
can see that a logarithm is also allowed. The variation of a logarithm is a constant, and in 
position space, this is just a contact term. In other words, we can allow c 7^ in the right 
hand side of (14.171) . Another way to see this is to consider the Fourier transform of the 
dilatation constraint. When we substitute x — idk we are implicitly integrating by parts. 
In general we neglect the surface terms because they are not singular. In the case we are 
considering, one can see that these terms are non-zero, hence c 7^ 0. 

In order to fix the combination of logarithms we can impose permutation symmetry as well 
as a good OPE expansion. The OPE expansion in position space says that {000) ~ -3-^ 



following solution is allowed 



as X23 — 0. This translates into (000) ~ -U^ as A;i — )■ We then find that only the 



f{k,,k2,k3)^log{k^ + k2 + ks) (4.18) 



Note that in this case the Fourier transform has dimension minus one, and then the special conformal 
generator is given by the simple expression in (|4.13p . 

This OPE requirement is imposed up to contact terms. Thus, for example, a term of the form logfc2, in 
this limit gives us a S'^{xi2) (since it is independent of fci) and is consistent with the OPE requirement, 
which is only imposed at separated points. In other words, when we expand (|4.18p for small fci we get 
2 log k2 + + • • • , we drop the first term and the second leads to the correct OPE. 
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It is possible to check that this is also the Fourier transform of the usual position space 



expression, — \ — 3-. It is also possible to show that one has simple solutions when operators 

^12^13^23 

of A = 2, 1 are involved. This is done as follows. After we obtain f l4.18p we can express 
the Fourier transform of the three point function as / = 117=1 kt'~^9- Then g has scaling 
dimension zero, and the special conformal generator on each particle acquires the form 

3 



b.k 



\i=i 



-(A-2)(A-1) d' 
\k[^ ^ d\k\ 
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(4.19) 



We then see that for A = 1, 2 the computation is the same as what we have done above. If 
all operators have A = 1, then the answer is = 1 or /a=i = k^^g^- When some operators 
have A = 1 and some A = 2 we cannot use permutation symmetry to select the solution, 
but it should be simple to find it. 



4.4. Constraints of special conformal invariance on conserved currents 

The constraints for special conformal invariance in momentum space are given by (14.121) . 
Here we would like to express the constraint of special conformal invariance in the spinor 
helicity variables. We would like to express the special conformal generator in terms of a 
simple operator such as (I4.13p . The operator we want to consider is the current in Fourier 
space, multiplied by a polarization vector proportional to = In fact, just multi- 

plying by this vector has a nice property, it leads to an operator of dimension minus one, 
since the Fourier transform of a conserved current has dimension minus one, and this choice 
of polarization vector does not modify the scaling dimension. If J is the conserved current, 
we take ^~.J and we act with the operator (I4.13p . The lambda derivatives can act on ^ 
and also on J, when they act on J, we can express them in terms of k derivatives. After a 
somewhat lengthy calculation, one can rearrange all terms so that we get the action of (I4.12p 
on the current, plus a term proportional to the divergence of J, or k.J. More explicitly, we 
find^s 

The first term in the right hand side vanishes due to the special conformal generator. The 
second term in the right hand side involves a longitudinal component of the current. One is 
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M\- - S) and D are defined in 
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tempted to set that to zero. However, we should recall that, inside a correlation function, 
we get contact terms at the positions of other charged operators. These terms are simply 
given by (the Fourier transform of) the Ward identity 

n 

k{{J\h)02{h) ■ ■ ■ On{K)) = -J2 Ql{02{k2) ■■■Ol{ki + h)--- Onikn)) (4.21) 

1=2 

Where Qi is the charge of the operator Oi. These are lower point functions. The conclusion 
is that there is a simple equation we can write down, by acting with the special conformal 
generator in spinor helicity variables, fl4.13p . 

Note that in position space, we normally impose the special conformal transformation 
at separated points. In other words, we do not consider local terms. A local term will 
contribute to the three point function in position space as 

{0{x)0{y)0{z)) Local ~ V[S%x - y)]f{x - z) + cyclic (4.22) 

V is an arbitrary differential operator, which, when integrated by parts, will just yield an 
analytic function of the ks (like kl, kfk^, ...). Upon Fourier transforming the first term, we 
see that its form will be [analytic piece] x F{k^) + cyclic, where F is the Fourier transform 
of /. So, any piece in the three point function that is analytic in two of its variables, like 
kl, kikl, k^/ki, corresponds to a local term. Something like kik2 is analytic in k^ but not 
in kl and k2, so it is non-local. 

Let us see how this works more explicitly. We can start with the two point function 

{r.J{ki)r-J{k2)) = S'iki + k2)^^ (4.23) 

Here the right hand side of the Ward identity vanishes and indeed, this function is annihilated 
by gH. 

Now we consider the three point functions of currents J" which are associated to a non- 
abelian symmetry. In the bulk they arise from a non-abelian gauge theory. The Ward 
identity is given by 

kl{Jnki)J!;ik2)J[ik,)) = r'%J,{ki + k2)Ji{k,)) - {J,{k2)Mk, + kl))] (4.24) 
Where the color factor was stripped off the two point correlators 

One unpleasant feature of this equation, (|4.24|) , is that the currents in the right hand side are evaluated at 
a shifted momentum, so it is not trivial to express this in terms of the A and A variables. In this particular 
case, this is not a problem since the two point functions can be explicitly computed, but this might lead 
to a more complicated story in the case of higher point functions. 
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Just as a check, let us compute the three point function for a gauge theory with a 
Yang Mills action in the bulk. Since the gauge field is conformally coupled, we can do the 
computation in flat space. We compute this correlator between three gauge fields in flat 
space, all set at t = and Fourier transformed in the spatial directions. We have the usual 
jabc non-abelian coupling in the bulk. In Feynman gauge, the final answer is 

{A%{h)A;i{k2)A%{h)) oc 6'{h + h + ^3) |^^||^^||^^| ;^r^^^(fef - ^'2') + cyclic] (4.25) 
where E = l^/l-^"^ 

We now multiply by ^~ for each particle to compute the correlator and we get 

{^^.A-^{h)^^.A'^^{k,)^^.A-^{h)) (X 6%h + h + ^3)r'^^'^^p;^||^|^ (4-26) 

Note that the expectation value we started from, f l4.25p . is not gauge invariant. On the 
other hand, once we put in the polarization vectors and we compute the transverse part, 
as in fl4.26p . we get a gauge invariant result (under linearized gauge transformations). Note 
now that this expectation value is related to the currents, via a formula similar to (14. 3p . 
which introduces extra factors of the two point function, which here are simply factors of 
\k\. Thus we find 

(ei-^"n^i)6.^"H^2)6-^°n^3)) ocfciA;2A;3(ei.A'^H^i)6-^"Hfc2)e3.^'^nfc3)) (4.27) 

Now let us check that this expectation value obeys the conformal Ward identity, with the 
operator (14.130 . The action of (I4.13p on the first current is 



... d' 



d\ldX\ 



(1,2)(2,3)(1,3) 
(1,1)(2,2)(3,3) 



A2aAla(2,3)(3,l) A3aAla(l,2)(2,3) 



1,1)2(2,2)(3,3) (1,1)2(2,2)(3,3) 
2Ai,Ai,(l,2)(2,3)(3,l)" 



(1,1)3(2, 2) (3, 3) J ^^-^^^ 
We now use the Schouten identity - a consequence of the fact that the spinors live in 
a 2D space - to simplify this further. Expressing A2a in terms of Aia and Aia we have 
(1, l)A2a = -(i, 2)Ai, + (1, 2)Ai, and thus: 

A2aAi, (2, 3) (3, 1) Ai„Ai, (1, 2) (2, 3) (3, 1) Ai„Ai^(l, 2) (2, 3) (3, 1) 



:i, 1)2(2, 2) (3, 3) (1, 1)3(2, 2) (3, 3) (1, 1)3(2, 2) (3, 3) 



(4.29) 



This is done as follows. The three point function in the bulk is the term in square brackets. We attach 
the propagators. We write the energy conservation condition as / dt'e^* ^ . Then we integrate over fc° 
to localize things at t = 0. We can deform the integration contour and we only pick the residues of each 
of the poles of the propagators which give rise to the factor in the denominator. We close the contour up 
or down depending on the sign of t' . In each case, integrating over t' gives us the factor of 
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We can do the same for X^a and then all that remains is a term proportional to Ai^Aia, given 
by 



laa 



(1,2)(2,3)(1,3) ' 
dXldXl L (1,1) (2, 2) (3, 3), 



(l,i)^(2,2)(3,l)[<^'2^<2'^^^^'^^-^^'^)<^'2^^2'^^] 

(4.30) 



Using the momentum conservation condition we can express "cross-products" of the form 
(m, n) for m 7^ n in terms of other brackets (the details are worked out in an appendix) 
through {ki + k2 + k^){m,n) = —2{m, o) (o, n) , where m ^ n ^ o and then the term in (I4.30p 
is 

(1,2)(2,3)(3,1) - (1,3)(1,2)(2,3) = -(2,3)2(A;2 - h) (4.31) 

Putting the pieces together, this is the expected contribution from the Ward identity 
3 



1=1 

=r"' l7%T-^2#^[(2,2) - (3,3)] + cycUc] (4-32) 
^ 1(1,1)2 (2,2)(3,3)^^ > / \ , /J y ( 



=r'^|^6-6[^2-A;3] + c2/c/^c| 



The expectation values that we have computed for gauge fields in de Sitter can also 
be computed in flat space, since gauge fields are conformal invariant So, we are simply 
computing correlation function of gauge invariant field strengths in but on a particular 
spatial slice. We are putting all operators at t = 0. We can think in momentum space 
and consider the operators Fahii — 0, ^), Pbh^ ~ 0' ^) where we Fourier transformed in the 
spatial coordinates but not in the time coordinate. Given k for each operator, we can define 
A and A via (13. 3p . We then can write the operators we considered above as: 

In both of these expressions, when we write Fab, or F^^^ we mean the self dual and anti-self 
dual parts, but the indices are summed over with the indices of the indicated A's. These 
expressions involve contractions that are not natural in flat space, but are reasonable once 



This is true for the tree correlators we are discussing but it is not true if loop corrections are taken into 
account. 
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we break the full Lorentz symmetry to the rotation group. These operators are set at t = 0, 
and in Fourier space in the spatial section, with momentum k}^ 

One can write higher derivative operators that give rise to three point functions that are 
annihilated by the special conformal generator fl4.13l) . These operators would be aTr[F^] 
and bTr[FF^]. Now it is important to put in the dS metric. The wavefunctions for A are 
still the same as those in flat space, if we compute the three point functions perturbatively. 
These give the following three point functions^" 

(J-+(l)/'+(2)J-+(3)) oc i2rrmY: + ^ ^)r" \^^^ j) '+ (2,'|+ (3! 3))3 ^^'^^^ 
(J"'-(l)j'''-(2)J^'-(3)) oc {2t:)H\Y ki){a-ih)r'' %, ^I^^'-^^^'/^-vn. (4-34) 



The result (14. 26 p . which comes from the usual Yang Mills term can be converted into a 
correlator of curents of the form 

(r'+(l)/'+(2)J^'+(3)) - '^'^f^^ 



M"'+(l)M^'+(2)M^'+(3) 



oc.-(.^, + 4 + W— l^llfif (4.35) 

is completely analytic in momentum space, and could be viewed as arising from a factor 
in the wavefunction of the form \1/ ~ ^Tr[A/\A/\A] ^ However, we would need a term in the 

wavefunction with the opposite sign to remove the correlator. Thus, although it 

looks like a local term, it does not seem possible to remove both the + + + and the 

correlator with the same factor. On the other hand, (14.331) is definitely non analytic in 
momentum, due to the singularity. 

The current correlators are derivatives of the wavefunction. The expectation values of A 
can be obtained from them. In that case the parity violating h term in (I4.33P drops out. 



Note that both the self dual and anti-self dual parts of F contribute to each of the helicities. The reason 
is that we have defined a four momentum (|/c|,A:) in order to define A, A. (We could also have reversed 
the sign of the zeroth component to |fc| — — |fc|, which exchanges A o A). With this definition, the time 
component of this four momentum is not necessarily equal to the total four momentum of on shell waves 
coming in or out of the F insertions (it can differ by a sign) . 

These couplings do not require a non-abelian theory. They are antisymmetric in the Lorentz indices, thus 
they require an antisymmetric tensor. Thus, we could have three abelian gauge fields F^ and then the 
cubic couplings eijLtr[F^ F^ F^], where the trace is over the Lorentz indices. 
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4.5. Constraints of special conformal invariance on the stress tensor 



In this section we consider the constraints of conformal invariance in momentum space 
for the stress tensor. 

We multiply the stress tensor by a convenient polarization tensor e~Tjj(A;), with e^- = 
In order to study the action of the special conformal generator it is convenient to 
define an operator containing an extra power of k as 

T = , e,. = = (4.36) 

The power of k was chosen so that the special conformal generator has the simple expression 
given by f l4.13p . The expression for is that one that would give a naturally normalized 
tensor, when we take the reality conditions into account. Again, this operator does not quite 
annihilate the correlator, but it produces a term involving the Ward identity in the right 
hand side. Although more laborious in terms of manipulations, the general ideas are the 
same as in the current case so we will be more brief in the details of the conformal symmetry 
check. 

We find that fl4.13p acts on the stress tensor as 

jIh = ^[4M^-, - S^^D]Tu - 3^b,e;^ki{Tij + T,i) (4.37) 

The first term is what we expect from fl4.12p for the stress tensor, and it vanishes. The 
second term can be computed by using the Ward identity. Again, such terms are analytic 
in some of the momenta. So if we disregard terms that are analytic in the momenta, we can 
drop also the term involving the Ward identity. 

Let us first ignore this Ward identity terms and compute homogeneous solutions of the 

equation. Let us consider first a general three point function of the stress tensor. 

Such a general three point function is given by 

{f-f-f-) = 5%J2 2) (2, 3) (3, 1)] V(^i, A:2, k,) (4.38) 

With / a symmetric function of dimension minus six. After some algebra, using Schouten 
identities, we get 

E(^T'^^3;^ 2) (2> 3) (3, 1)) V(fci, A:3)] = -2|(1, 2) (2, 3m 



h.OT- = b.O 



+ K1,2)(2,3)(3,1)] 



k\ + cychc (4.39) 
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Although the are hnearly independent, the ks are not. A convenient way to rewrite fl4.39p 
is to choose special conformal transformation parameters 6* to project out a few components. 
Let us take 6* ~ (-^2-^3 + -^3-^2)*) example. This combination was chosen so that the time 
component of b is zero. We find 

{X2X3).d[i )/] = (1, 2)^(2, 3)^(3, 1)^ {4(9,, - 9,3)/ + hidl - dl)f - k,{dl - dl)f} 

^ = 4(9., - d,,)f + k,{dl - dl)f - k,{dl - dl)f (4.40) 

It is straightforward to check that the gravity result we obtained from a in interaction 
contribution gets annihilated by this operator. Such a contribution is simply / = (^1 + ^2 + 
^3)"^. The Einstein contribution is not annihilated. It gives a nonzero answer that matches 
the expected answer from the Ward identity. One can solve the equation fl4.40p and its two 
other cyclic cousins by brute force. One finds the expected solution, mentioned above, plus 
a new solution that has the form 



' [(A;i + k^- k^){k^ + ks- k,){k3 + k,- k^W ' ■ ' 

This new solution does not have the right limit when ki — )■ 0. Namely, if we start with 
Tij{ki), when the momentum goes to zero we do not expect any singular term when — )■ 0. 
In fact, an insertion of Tij{k = 0) corresponds to a constant metric or a change of coordinates. 
So, in fact this limit has a precise form. On the other hand if we look at this limit in fl4.4ip . 
we find that / ~ l/fcf, which is too singular compared to the expected behavior. 

Let us now turn our attention to the h correlator. We first write an ansatz of the 

form 

(f -T-f +) = 6'iJ2 2) (2, 3) (1, 3)]'gik,, k^, k,) (4.42) 

where 5^ is a homogeneous function of degree six We can now use a trick to get the 
homogeneous solutions of the special conformal generator. We note that if we exchange 
A ^ A, then the sign of \k\ is changed, but k does not change. Now, the ansatz for f l4.42p 
differs from f l4.38p precisely by such a change in the third particle. Thus, the two solutions 
for g in f l4.42p are simply given by the two solutions for / but with k^ — )■ —k^. More 

In analogy to flat space, one is tempted to write this in terms of [ (2^3x1 3) '^^^^ easy to do using the 
identities in appendix |B1 but we found simpler expressions in terms of (|4.42p . 
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explicitly, the two solutions are 

1 

^ =Pi + h + h){h + h- + h- hW 

However, now both solutions are inconsistent with the small limit. The first solution has 
a problem when fci — )■ and the second when k^ — t- 0. Thus, both are discarded since these 
limits are too singular. Even though this trick of exchanging A ^ A was useful for generating 
solutions of the homogeneous equation, the full results for the correlators are not given by 
such a simple exchange. 

In conclusion, we have shown that there are no other solution of the conformal Ward 
identities beyond the ones we have already considered. It remains to be shown that the 
Einstein gravity answer obeys the conformal Ward identity. One can check that the Einstein 
gravity answer is not annihilated by the operator fl4.13p . It gives a nonzero term. This term 
is indeed what is expected from the Ward identity for the stress tensor, which is the second 
term in fl4.37p . Of course, this is expected since Einstein gravity has these symmetries. The 
relevant expressions are left to appendix lEl 

5. REMARKS ON FIELD THEORY CORRELATORS 

In this section we compute the free field theory three point correlation function for scalars 



17j in position space. Here we work 



and fermions. This is very similar to what was done in 
in momentum space. We will compare these expressions to the gravity ones computed above. 
The idea is that by considering a theory of a free scalar and a theory with a free fermion 
we obtain two independent shapes for the three point function of the stress tensor. Since 
we have computed the most general shapes above, this will serve as a check of our previous 
arguments. In addition, the momentum space expressions for the correlators might be useful 
for further studies. The two correlation functions that we obtain for scalars or fermions are 
parity conserving. Our results indicate that there can be field theories that give rise to the 
parity breaking contribution. Such field theories are not free, and it would be interesting to 
find the field theories that produce such correlators We will concentrate here on the free 



There are bounds similar to the ones derived in [32i] for the parity breaking and parity conserving 



coefficients that appear in a three point function. The parity conserving bounds were considered in 38 1. 
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theory case^^. 

The computation is in principle straightforward, one simple has to compute a one loop 
diagram with three stress tensor insertions. The only minor complication is the prolifer- 
ation of indices. To compute the diagram itself one needs to use the standard Feynman 
parametrization of the loop integral. In particular, it is convenient to use the following 
Feynman parametrization: 

1 2 da d/3 

ABC ~ Jo {A + aB + ^Cf ^ ' ' 

The final expression will have several contractions of polarization tensors with 3-momenta, 
e.g. eije|fce|j, ^^ik'^^ljkfef^k'^el^^, etc. Then one needs to use the expressions presented in the 
appendix to convert these to spinor brackets. The final answers have a simpler form than 
the ones with the polarization tensors. 

We treat first the scalar case and then the fermion case. 

5.1. Three point correlators for a free scalar 

The stress-energy tensor for a real, canonically normalized scalar field is 

T,,{x) = h,<f){x)d,<j){x) - ]mdj(l){x) - l6,,d^<f)\x) (5.2) 

4 4 o 

The two point function is, up to the delta function: 

iT*T-U = ^ (5.3) 



These bounds can be derived by considering a thought experiment where we insert a stress tensor at the 
origin with some energy and then we look at the angular dependence of the energy one point function, 
as measured by calorimeters placed at infinity. The stress tensor at the origin has spin ±2 under the 
S0(2) rotation group of the spatial plane. The energy one point function as a function of the angle is 
{£{0}) ^ e*''^ + 1 + e"***^, with coefficients that depend on the parts of the stress tensor three point 
functions that we have characterized as coming from (W^)^, R, (W^)^. 
^■^ Similar calculations involving the trace of the stress tensor were considered in [22| . 
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The three point functions are 

{T+{h)T+{k2)T+{h)) 
(ki + k2 + k3f 



kl + ki + kl (/ci/c2/^3) 



+ 



128 



{ki + k2 + ks) 



64 ■ 2{ki + k2 + k3)^ 

{kik2 + kiks + k2kz) kik2ks 



(ki + k2 + ks) (ki + k2 + k^Y 



{T+{k,)T+{k2)T-{ks))^ 
(ki + k2- k^f 



2-3 _|_ jL3 _|_ iL3 

nj-^ \^ rv'2 \^ 'i'3 



64 



128 



(/ci + k2 + ka) - 



{kik2 + kiks + k2k3) 



kik2k3 



{ki + k2 + ks) {ki + k2 + kaY 



((1,2)(2,3)(3,1))^ 
(5.4) 



((i,2)(2,3)(l,3)r 

'^l'^2'^3 

(5.5) 



5.2. Three point correlators for a free spinor 



The stress tensor for a complex Dirac spinor is given by 



The two point function is 



128 



The three point function is given by: 

i-S I 1,3 I jL3 



(T+(A;0T+(A;2)T+(A;3))^ 

{ki + /c2 + h? 



{kik2k3f 



64 



{ki + k2 + ks 



64 



(/ci + k2 + ks) - 



{kik2 + kiks + k2k3 



kik2k3 



{ki + k2 + ka) {ki + ^2 + h)' 



{T+iki)T+ik2)T-iks))^ = 
{ki + k2- k^f 



_L _L h'^ 



64 



64 



{ki + k2 + ks) - 



{kik2 + kiks + k2k3 



kik2k3 



(/ci + k2 + k^) {ki + k2 + k^Y 



(5.6) 



(5.7) 



((l,2)(2,3)(3,i))^ 

J,2l2l2 
"'1'^2"'3 

(5.8) 



((l,2)(2,3)(i,3))^ 

h.2h.2h.2 

"'1"'2"'3 

(5.9) 



5.3. Comparison with the gravity computation 

We sec that these results contain the general shapes discussed in gravity, but they also 
have an extra term proportional to ^ /cf . This is a contact term. Namely, it is non-zero only 
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when some operators are on top of each other. In position space we get a delta function of 
the relative displacement between two of the insertions. These terms are easily recognized 
in momentum space because they are analytic in two of the momenta. These contact terms 
represent an ambiguity in the definition of the stress tensor. There is no ambiguity in 
taking the first derivative with respect to the metric. However, these contact terms involve 
a second derivative with respect to the metric. So if we define the metric as g = e"^ and we 
take derivatives with respect to 7 we are going to get one answer. If we took (7 = 1 + 7' and 
took derivatives with respect to 7' we would get a different answer. In fact, we have the 
same ambiguity in the gravity results if we define 7^'^- = •jij + ^'jwyij. In that case the two 
results will differ precisely by such a term. It is interesting to note that the non-gaussian 
consistency condition discussed in {sj] does depend on this precise definition of the metric, 
since a constant 7 gives rise to different coordinate transformations depending on how we 
defined 7. The one derived in js] holds when the metric is defined in terms oi g = . 

We can note that if we have scalars and dirac fermions, then we have to sum the 
two contributions to the three point functions that we have written above . li = 2n^ , 
then we see that the term going like 1/ cancels. This is the contribution that comes from 
a term in the bulk. This combination is also the one that appears in a supersymmetric 
theory. In fact, in a supersymmetric theory the three point function of the stress tensor does 
not have any free parameters'^. It is the same as the one given by the pure Einstein theory 
in the bulk, which does not contain any 1/ terms. This is related to the fact that in four 



flat dimensions supersymmetry forces the + + + and amplitudes to vanish [2 



6. DISCUSSION 



In this paper we have computed the possible shapes for non-gaussianity for gravitational 
waves in the de Sitter approximation. Though three possible shapes are allowed by the 
isometries, only two arise in de-Sitter expectation values. The parity violating shape con- 
tributes with a pure phase to the wavefunction and it drops out from expectation values. 
The two parity conserving shapes were given in equations (12.111) . (12.181) . One of these shapes 
is given by the Einstein theory. The other shape arise from higher derivative terms. Under 

Of course, to have supersymmetry we need to consider an AdS, rather than a dS bulk. 
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general principles the other contribution can be as big as the Einstein term contribution. Of 
course, in such a case the derivative expansion is breaking down. However, the symmetries 
allow us to compute the three point function despite this breakdown. This is expected for 
an inflationary scenario where the string scale is close to the Hubble scale. This requires a 
weak string coupling, so that we get a small value of H'^/Mpi ~ g'^. One of these shapes is 
parity breaking. These three point functions of gravitational waves are expected to be small, 
having an fNL-gravUy = ^^^^p of order one. In a more realistic inflationary scenario, which 
includes a slow rolling scalar field, then we expect that these results give the answer to lead- 
ing order in the slow roll expansion. It would be interesting to classify the general leading 
corrections to the graviton three point function in a general inflationary scenario. This can 



probably be done using the methods of 



-[6|. Here by assuming exact de Sitter symmetry 



we have managed to compute the correction to all orders in the derivative expansion. 

We have presented the result in terms of the three point function for circularly polarized 
gravitational waves. We used a convenient spinor helicity description of the kinematics. 
These spinor helicity formulas are somewhat similar to the ones describing flat space am- 
plitudes. It would be interesting to see if this formalism helps in computing higher order 
amplitudes in de Sitter space. The problem of computing gravitational wave correlators 
in de Sitter is intimately related with the corresponding problem in Anti-de Sitter. (The 
two are formally related by taking R'^g — )■ —R\4si where R are the corresponding radii of 
curvature). Thus, all that we have discussed here also applies to the AdS situation. The dS 
wavefunction is related to the AdS partition function. In this case all three shapes can arise. 
The parity violating shape arises from the j W'^W term in the action. These three point 
correlators, for Einstein gravity in AdS, were computed in [t]. It would be interesting to 
see if the spinor helicity formalism is useful for computing higher point tree level correlation 
unctions. It is likely to be useful if one uses an on shell method like the one proposed in 
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4l| In the spinor helicity formalism that we have introduced, we have defined the 
"time" component of the momentum to be \k\. The choice of sign here was somewhat arbi- 
trary. When we are in four dimensional flat space, there is a simple physical interpretation 
for the results we get by analytically continuing to —\k\, as exchanging an incoming into 



Raju's proposal [40|, |41[ is only for > 4 dimensions, but it might be possible that something similar 
exists in = 4. 
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an outgoing particle. It would be interesting to understand better the interpretation for 
this analytic continuation of the correlators we have been discussing. This continuation was 
important in 
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41|. 



These computations of three point functions in de Sitter or anti-de Sitter are intimately 
related to the computation of stress tensor correlators in a three dimensional field theory. 
In fact, the symmetries are the same in both cases. Therefore the constraints of conformal 
symmetry are the same. The physical requirements are also very similar. The only minor 
difference is whether we require the two point function to be positive or not, etc. But 
in terms of possible shapes that are allowed the discussion is identical. Thus, our results 
can also be viewed as giving the three point correlation functions for a three dimensional 
ield theory. The position space version of these three point functions was discussed in 
3, I42I44I. For some three point functions the position space version is much simpler. 



On the other hand, for the stress tensor, the position space correlator has many terms 
due to the different ways of contracting the indices jl7]. The momentum space versions 
we have written here are definitely shorter than the position space ones. They are a bit 
convention dependent due to the contact terms. Thus, they depend on precisely how we are 
defining the metric to non-linear orders. Here we have made a definite choice. An elegant 
and simple way to write correlation functions is to go to the embedding space formalism 
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47|. It is likely that one can obtain relatively simple expressions for the three point 



correlators using that formalism. On the other hand, the momentum space formalism might 
be useful for constructing conformal blocks, since, in momentum space, there is only one 
state propagating in the intermediate channel. 
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Appendix A: Expression for the three point function in terms of an exphcit choice 
for polarization tensors 

As we are studying three point functions, there is a way to define polarization tensors 
that are similar to the usual "x" and "+" of General relativity. We call them X and P 
here, so as not to cause confusion with the helicity labels + and — . The choice of helicity 
states is based on the little group of an Euclidean 3D CFT. Basically, one takes two possible 
polarizations, P and X, as functions of a vector orthogonal to the plane of the triangle and 
of a vector orthogonal to one of the momenta we are looking at. So, using the notation 
defined in figure [H we have that 

e^P = 2{z,z,-uTuJ) (Al) 
4'- = 2{uTz,+z.uJ) (A.2) 

And the previously discussed + and — polarizations will be given by ± = P ± iX. P and X 
are the polarizations known as + and x in general relativity, but we choose to use different 
labels so that the former is not interpreted as positive helicity by mistake. 




FIG. 1: The 3-momenta and the auxiliary vectors used to define the polarizations. 

We list here the results for the non-gaussianities due to the Einstein term and the Weyl 
term. The relevant pieces are labeled by the polarization choices PPP and XXP. We 
always take particle three to have polarization P. The other structures are obtained by 
cyclic permutation. There are no PPX and XXX structures because they break parity. 
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since z flips under parity so that X is odd and P is even. We use here the notation 
A;2, k^) = 2{klkj + kfkl + k^kj) - {kf + ki + kj). 



H 



-1 



ki + k2 + ks 



kik2 + kikz + k2kz 



MpiJ A{kik2ks) 
kik2k3 



j{kr,k2,ks){ykt+6yk^k 



j=i 



i<j 



ki + k2 + ks {ki + k2 + k^Y 



(A.3) 



{llllli)R={'^^f5'[Y.k 
ki + k2 + k3- 



H 



1 



Mpi) {kik2k3)^ 



J{ki,k2,k3) 



ki + k2 + k3 {ki + k2 + k^Y 



(A.4) 



Pi 



2160 



a; (A;i + A:2 + A;3)4(A;iA;2A;3)' 



;J{ki,k2,k3) = 



/J V 270(A;i + A;2 - A;3)(A;2 + k^ - ki){k3 + ki - k2) 



(A;i + A;2 + A;3)3(A;iA;2A;3)2 



(A.5) 



(7fei7fe2 7fe3)M^3 = -(7fei7fe2 7fe3)i^^ 



(A.6) 



-l) 










i-i oJ 


[l J 







Appendix B: Details on the spinor helicity formalism 

Here we summarize some conventions that we have used. 

• Metric: 77^^ = diag(-l, +1, +1, +1); 

• Sigma matrices: a^"'b — {—5"'b,a^°'b); 

• Scalar product: p.q = —2{Xp,Xq)(Xp,Xq); Energies: = p= — (Ap, Ap) — — CafcApAp 

• Polarization vectors used for the expressions in the appendix (normalization is not 

A'^A^ 

the same as the one used in the paper for the stress tensor): — ; — and 



k 



i+aa 
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— * — * 

Starting from a three momentum k, we define a four momentum k*^ — {\k\,k). This 
obeys k^k^ — 0. This defintion can be done for dS, AdS, or a three dimensional CFT. Note 
that k^ = a^ab^''X"', where a* are the Pauh matrices with an index lowered by eac- These 
matrices are symmetric. Thus we conclude that if we exchange A'' -H- A'' we keep the value 
of k, but we change the sign of the energy = k. This change is not consistent with the 
reality conditions which are (A*)" = ^ha}"^, (A*)^ = ~^ha^"'- However, if we just forget about 
the reality condition, then the exchange of A ■<-> A is allowed. Note also that the following 

— * 

symmetry is consistent with the reality conditions and it reverses the sign of k but does not 
change the sign of k. Namely, the exchange A' = A, A' = —A. This is useful for the two 

point function, where the momentum conservation condition forces the spatial momenta to 
be opposite. Note that in some formulas we write /cj, which means |A;|. This can also be 
written in terms of ki — —(A, A). 

— * — * 

For a given 3- momentum k — {ki,k2, k^) if we define |A;| = kg then one can for example 
take the following explicit choice of spinors, assuming that the reality condition is satisfied: 

Let us now summarize a few identities that are useful for the treatment of the three point 
function. For the case of the three point function the 3-momentum conservation condition 
reads 

F ^ 
XfXt + A»A« + XfXi = --e- , E = k, + k2 + ks^- ^(A„, A„) (B.2) 

n=l 

Where the coefficient two is determined from contracting the a and d. We can contract this 
expression with, say, Ai and A2. The purpose of that is to derive an expression for an object 
of the form (m, n), which has no interpretation as an energy, ii m ^ n. By doing that we 
find that 2(1, 3) (3, 2) = E{1, 2) and hence, we can write the general expression: 

_ ^ 2{m,o){d,n) 
{m,n) = - 2^ (B.3) 

Where m ^ n ^ o. Also, note that the 4-momentum product of two distinct 4-momenta, 
in terms of the three energies and the total energy, is given by: 

- 2{m,n){fh,n) = k'^k''^' = + k^.k" = ^{ko - k^ - kn)E (B.4) 
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And we also use the Schouten identity, which is useful to write a given spinor A in terms 
of two reference spinors and ^ 

(e,^)A'^=(e,AK-(/i,A)r (B.5) 

Whenever we have an expression in terms of angle brackets we can write it completely in 
terms of (m,n) brackets and the km = —{m,rh) brackets. Useful identities to do so are (for 

m ^ n) 

I -\ o) \ko ~\~ kn — km) / _ _\ [ko ~ km — kn)E , , 

{m,n) = - , (m,n) = ^ (B.6) 

{n,o) 2 A{m,n) 

1. Some expressions involving polarization vectors 

Now, let us calculate all the possible contractions of polarization vectors and momenta 
of different helicities: 

(B.7) 
(B.8) 
(B.9) 



Sm Sn 


~ h h 
<^m<^n 




~ h h 
'^m'^n 


Sm Sn 


h h 



As for contractions of momenta with polarization vectors, we have 
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k„.C = - '<'"■"><'"■"> (B.IO) 



2(m,n) (171,12) ,^ 

These can also be used to convert the gravity expressions into expressions in the spinor 
helicity variables since the gravity polarization tensor is e^j ~ C,iC,j- 



Appendix C: Comments on the parity breaking piece of the two point function 



As pointed out in 26| , the gravitational wave two point correlation function (or grav- 



itational wave spectrum) can be different for the two circularly polarized waves without 



As there is no time component of the polarization vector, k°^°^€aa ~ k^e^^ — kiCi so we are only taking the 
space components into account. 
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breaking rotation symmetry. In fact, a bulk coupling of the form J f{(f))WW is enough to 
produce this. This mechanism requires an inflaton. One can ask whether a parity breaking 



28j. Here 



two point function is possible in de Sitter space, as some authors have suggested 
we make some comments on the parity breaking pieces of the two point function of the stress 
tensor. 

The summary is that parity breaking terms are allowed in the gaussian part of the 
wavefunction of the universe, or in the two point function of CFT's. However, such terms 
are local, and contribute with a phase to the wavefunction. Thus they do not lead to different 
amplitudes for left and right circular polarizations. 

Let us start by discussing this from the wavefunction of the universe point of view. From 
that point of view the question is whether there can be a parity breaking two point function 
for the stress tensor. One is tempted to say that the answer is no. The argument is the 
following. The stress tensor is in a single representation of the conformal group, thus its two 
point function should be uniquely fixed. In fact, this is correct if we consider the two point 
function at different spatial points. However, there can be a parity breaking contact term. 
In order to understand this, let us discuss first the case of a current, or a gauge field in the 
bulk, and then discuss the case of the graviton or the stress tensor. 



1. Parity breaking terms in the two point functions for currents 



These were discussed in the AdS context in [48|. We just summarize the discussion here. 



The Fourier transform of the conserved current two point function is (in a 3D CFT) 

{Ji{h)Jj{k2)) = 5%h + h) - hkj)\k\-^ + ^e,,,fci,i] (C.l) 

This is consistent with conformal symmetry. It is annihilated by f l4.12p . The 9 term breaks 
parity. Since this term is analytic in the momentum, it gives rise to a contact term in position 
space, a term proportional to ieijid^iS^ {x — y) . If we couple the current to an external source 
and compute "^[A] = Z[A.J], then we are just adding a local term to the wavefunction 
of the Chern-Simons form ipeiA) = e^^ '^^^'^'^^ipg=Q{A) This is what we would get in 
a dS situation if we have an ordinary 6 term in the bulk. In other words, if we have a 
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For the non abelian case, we can complete this quadratic action into the full cubic one. 



42 



gauge field in the bulk with an interaction 9 J Tr[F A F], then the wavefunction contains a 
term proportional to the Chern-Simons action on the spatial slice. In a unitary (and gauge 
invariant) bulk theory this term has a real value of 6, so that it contributes as a phase in 
the wavefunction. Thus, when we compute the square of the wavefunction, this terms drops 
out. More explicitly, we can now compute the wavefunction in momentum space 

^{A) = exp{-A{k)A,{-k){6,,\k\ - ^ + Oe,^sks)} (C.2) 

We see that the 6 term is imaginary if we take 6 to be real. (We use that A{k)* = A{—k)). 
Then, if we compute [^/^(A)^ we find that the 6 term drops out if 6 is real. 

Now we can ask, could it be that some unknown unitary Hamiltonian produces a wave- 
function that contains a Chern-Simons part with a purely imaginary 67 The arguments 
leading to the Chern-Simons term in the wavefunction were based purely on demanding 
conformal symmetry and did not rely on any assumptions about the bulk Hamiltonian, or 
even its existence. All we are assuming is that we have a wavefunction that is conformal 
invariant. In particular, purely from conformal symmetry, the 6 term could be imaginary. 
An imaginary 6 leads to different amplitudes for the two circular polarization states of the 
gauge field^^. 

One problem with this is that the resulting probability amplitude is now not invariant 
under large gauge transformations. This is due to the fact that the Chern-Simons action 
shifts by a certain real factor under a large gauge transformation. Thus the wavefunctions 
produced in this way are not gauge invariant This argument is most clear in a non-abelian 
situation. 

However, if we ignore this problem, then we should also point out another issue. A 
different amplitude for left and right circularly polarized waves violates CPT invariance^°. 
This is most clear if we think about the observer in static patch coordinates. This observers 
sees de Sitter as static. For this observer CPT is a symmetry, it is not spontaneously broken 
by the background. However, CPT transforms -|- circular polarization into — . Thus these 

We are not assuming that we have an ordinary term in the bulk, but simply that some unknown dynamics 
gives rise to a 9 term in the wavefunction as in (jC.2p . 
This argument was suggested to us by E. Witten. 
^'^ It is not known whether CPT invariance holds in quantum gravity. In theories that have a dual CFT 
description, like the ones in AdS, CPT is a good symmetry, so it seems reasonable to assume that CPT 
will still be a symmetry of dS. 
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amplitudes cannot be different. Note that the wavefunctions that we discussed are the late 
times ones, the wavefunctions for fluctuations outside the horizon. On the other hand, the 
static patch observer probes the wavefunction inside the horizon. So, here we have assumed, 
by continuity, that if we get a parity breaking effect outside the horizon, then we should also 
see some effect inside the horizon. 



2. Parity breaking two point functions for the stress tensor 

Similar arguments can be used for the two point function of the stress tensor in momentum 
space. The only term we can write that breaks parity, by power counting, is 

{Tij{k)Tmn{-k)) odd ~ [{^imlh^jnk'^ + {l ^ j)) + {ui ^ u)] (C.3) 

This is a function of fcj and fc^, hence, it is analytic and corresponds to a local term in 
position space. 

In gravity, an analog of the 9 term is the topological invariant 

jTr[RAR] = j e'^'f"' Rf^Rf^ = j t^"'"'' R^'^^W^ = j WW (C.4) 

The last equality follows from the symmetries of the Riemann curvature. Adding this term 
as ^ J WW to the bulk action we get a contribution to the wavefunction of the form 

where u is the spin connection. It can also be written in terms of the Christoffel connection 
Let us check that this term indeed produces (IC.Sp . We expand the Chern-Simons term 
to quadratic order and obtain 

Scs ~ I e'^^Tldij:% (C.6) 



Using the first order expressions of the connection we find 



Scs (^ijl{drlsi- dslri)dldrlsj ^ (^ijlk'^l{k) sihli-k) si (C.7) 



■^^ In fact, we can use the relation between the two connections that comes from demanding that -D^e" — 0, 
which is Lo"^ — e^F^^i?"'' — df^e'^E"'' In this form it has the form of a GL{N) gauge transformation, 
ujfi = gV ^g~^ — d^gg~^ with g = ej^. Of course, gauge transformations are a symmetry of the CS action. 
Thus we get the same action in terms of both connections. The first (upper) and last indices of F are 
viewed as the "internal" GL(N) indices of the connection. 
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where we used that 7 is transverse, ks'jsi = 0. This indeed reproduces (IC.Sp . 

As in the gauge field case, if 6 is real, this term disappears from I^E'p. On the other hand, 
if 6 is imaginary, we do get an extra contribution to which leads to a different amplitude 
for left and right circularly polarized gravitational waves, as pointed out in Q]. 

Note that Scs in (IC.Sp is not invariant under large gauge transformations of the local 
Lorentz indices of the spin connection. 

Again CPT invariance forbids a different amplitude for left and right circular polarization. 

Note that all the remarks in this section apply to the case of pure de Sitter. In the case 
that we have an inflationary background, time reversal symmetry is broken by the inflaton 



and we can certainly have a parity violating two point function [26 1. 



Appendix D: Commuting through the delta function 

In this appendix we show that the action of the special conformal generator (14.121) on a 
correlator or expectation value of the form 6{P)M is equal to the action of the operator on 
A4. In other words, fl4.12p commutes with the momentum conserving delta function. The 
point is to understand how to get all the momentum derivatives through the momentum 
conserving 6 function. From now on all derivatives will be k derivatives. 

The delta function depends only on the sum of all momenta, let us call that P. Then the 
sum over all particles of Da where a runs over particle number has the form 

J2 Da5\P) ) M = Q[{b.dp)6'^]M (D.l) 

\ a / 

In order to derive this we have done the following. In each term the derivatives are with 
respect to ka, which end up dp when acting on the 6 function. Then the ka in D all sum 
up to P. Thus we have a term of the form P{dpdpS)A4. We then integrate by parts the 
derivatives to act on Ai in such a way that we get terms of the form in (ID.ip and also terms 
of the form P6{P)dpdpAi. Such terms vanish. Thus fID.ip is the total contribution from 
terms with two derivatives on the delta function. 

We can now consider terms which have only one derivative on the delta functions. There 
are terms coming from D. Let us consider those flrst. The flrst term in D contributes with 

2[9p/](26\fc^)9,. (D.2) 
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The second term gives 

2[dp^5] {{-2ki)(h.d,^) - 2V{k.du)] (D.3) 

The (lD.2p and the first term in (1D.3P give the action of the rotation generators on the term 
multiplying the 5 function. These would make the correlator vanish if it was rotational 
invariant. On the other hand, we have indices, thus the correlator is rotational covariant. 
The action of the rotation generators has to results in some action on the indices. These 
must arise from the action of M on the 5 function, producing the spin generators. Finally 
we have the last term, which adds up to the dilatation generator. Such terms combine with 
a derivative of a 5 function in f lD.ip and also a derivative that comes from the first term in 
f l4.12p . They altogether sum up to 

- (6.9p5)2 I ^(A, - 3) +'i-Y,^a.dA (D.4) 

\ a a ) 

The last term is computing the overall dimension of the term that multiplies the 5 function, 
the +3 is taking into account the dimensions of the 5 function. And the first term is the total 
dimension of the (Fourier transform) of the external states. Thus, if the answer is dilatation 
covariant these terms will vanish. The total dilatation eigenvalue kdk of the coefficient of 
the delta function is then 3 + X]a(^a " 3). 



Appendix E: Checking the Ward identities for the stress tensor 

The Ward identities come from the statement that the wavefunction of the universe is 
reparametrization invariant ^E'[5'jj] = ^[dij + V(ifj)]. This then implies that 







(E.l) 



The stress tensor correlators are defined by taking multiple derivatives of the wavefunction 
and then setting g to the fiat metric. By taking multiple derivatives of f lE.ip we get the 
Ward identity which looks like 

r 

di{Tij{x)Ti^^raAyi) ■ ■ ■Tl„,mn{yr)) = ^ /^J^^^^^ (x - y,) (T^^ (?/l) ■ " ■T;/,„/(y,) ■ • • T,„,„„ (y^)) 

s=l 

(E.2) 

where D is a first order derivative (acting on x or y^) and the indices are contracted in 
some way. These terms come from acting with the metric derivatives on the explicit metric 
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dependence in the covariant derivative, etc. Notice that if all the points in the left hand side 
are different from each other, then the right hand side is zero. In this section we will assume 
that all the points yi, ■ ■ ■ ,yr are different from each other. The precise form of the contact 
terms in the right hand side depends on the precise definition of the "derivative with respect 
to the metric" . If we define the stress tensor as derivatives of the from -Ar , then the Ward 



0. 



identities can be found in However, in our case, we defined the stress tensor correlators 
as derivatives of the wavefunction with respect to jij, where we write the metric as g = e'^, 
see (14.11) . This leads to slightly different expressions. The difference is only present as extra 
contact terms that arise when we use the chain rule T°"^® = -ttt = ^ 



One can keep track of these extra terms and write the precise version of the Ward identity. 
For the case of the three point function, after going to Fourier space, we get the simple 
expression 

{[henhm2^,T{hms^3nh)]) = 6.^3 {^m2-^3 + ^2mu3 + ^3Mi-^2} [2ki - 2ki] 

(E.3) 

Which, for is given by: 

{[heT{k,)]U2nk2)]U3T{h)]) = 8(1, 2) (2, 3)^(3, 1) ^^^ + ^(fc^fc^ff " ^^'^^ 

From the point of view of the operator f l4.13p . this expression (]E.4p should be multiplying 
a term that goes like {ki^)^^{...). The coefficient of the three point function is then fixed 
by comparing this with the result of acting on the three point functions of Einstein gravity 
with (Km . 

For a general three point function, given by ((1, 2) (2, 3) (3, l))^/(A;i, /c2, /^s), the 

action of KT^ is 



(^T'OX^^ [((1, 2) (2, 3) (3, l)yf{ku k,, k,)] 



[(1, 2)(2, 3)(3, 1)]2 (6 + 2kTd,^)d,.f - k\d,.d,.J 



+ 



+ .e^.-'^iM [2(1, 2) (2, 3)2(3, 1) [(2, 1) (3, 1) - (3, 1) (1, 2)]]^, J (E.5) 



32 This gets simplified thanks to the fact that {Tij{-k)[£,£,T{k)]) = 2f ^^fc^. 

■^■^ When there is no index in the derivative, it is understood that we are taking the derivative with respect 
to the energy, or \k\. We hope this does not cause confusion in the expressions here. 
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In order to derive this, we used Schouten to express A^, each as a function of A^ and 
A^. As is a complex vector, one can show that ie^-'^'^i^ = Ci- Then the third hne is 
proportional to The next step is to write the first piece as a function of ki and not of 
its components. The final result is: 
... 



dX-dXt 



[((1, 2)(2, 3)(3, l)ff{h, h)] = -2^1(1, 2)(2, 3)3(3, - ^2]^, J+ 



+ [(l,2)(2,3)(3,l)f 



k\ + cyclic 



(E.6) 



Although the ^s are linearly independent, the fcs are not. A convenient way to rewrite flE.6p 
is to choose special conformal transformation parameters If to project out a few components. 
Let us take 6* ~ (A2A3 + A3A2)* for example. The constraint is 

{X,X,).d[{ )/] = (1, 2)2(2, 3)3(3, If {4(4. - 43)/ + Hdl - dljf - k,{dl - dl)f] 

(E.7) 

For the h correlator, the derivation is similar. 
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